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SUMMARY 


Experience during the past 6 years of operation of the Langley 0.3-Meter 
Transonic Cryogenic Tunnel has shovm that there are problems associated with 
the operation and control of such tunnels. One problem has been an inability 
to provide long-term accurate control of test parameters. Additionally, the 
time required for transition from one test condition to another has proven to 
be excessively long. Due to the high degree of process cross-coupling between 
the independent control variables, man-in-the-loop operation has proven to be 
much less efficient than desirable in terms of liquid-nitrogen and electrical- 
power usage. For these reasons, studies have been undertaken at the Langley 
Research Center to model the cryogenic-wind-tunnel process, to validate the 
model by the use of experimental data, and to construct an interactive tunnel 
simulator with the validated model. Additionally, this model has been used to 
design closed-loop feedback control laws for regulation of the temperature and 
pressure parameters. 

The global mathematical model that has been developed consists of coupled, 
nonlinear differential governing equations based on an energy-state concept of 
the physical cryogenic phenomena. Although the fundamental process is quite 
stable, the tendency for continuous drifting from balanced equilibrium energy 
conditions is prevalent. The process equations, the simulation responses, and 
the experimental data are provided here for future reference. Also included 
are the control laws and simulator responses using these closed-loop feedback 
schemes. 


INTRODUCTION 

Historically, capital and operating costs have tended to keep transonic 
wind tunnels small, and the many problems encountered at high pressures have 
tended to keep operating pressures low. The net result has been that existing 
(ambient temperature) tunnels operate at Reynolds numbers which are far too low 
to insure proper simulation of the flow experienced in flight, particularly with 
regard to shock — boundary- layer interactions encountered on modern high-subsonic 
and transonic aircraft. 

Of the various ways of increasing Reynolds number that have been tried or 
proposed for transonic tunnels, cooling the test gas to cryogenic temperatures 
(150 K or lower) appears to be the best solution in terms of model, balance, 
and sting loads, as well as capital and operating costs (ref. 1 ) . Personnel 
at the NASA Langley Research Center (LaRC) have been stuc^ing the application 
of the cryogenic-tunnel concept to various types of high Reynolds number tran- 
sonic tunnels since the autumn of 1971 and, through extensive theoretical and 
experimental studies, have successfully demonstrated both the validity and 
practicality of the concept. As a result of this work, the first major tran- 
sonic tunnel especially designed to take full advantage of cryogenic operation 
is now under construction at LaRC. This tunnel, the National Transonic 



Facility (NTP ) , will provide an order of magnitude increase in Reynolds number 
capability over existing tunnels in the United States when it becomes opera- 
tional in 1982. 

It can be shown that for equal test Reynolds numbers and for any arbitrary 
maximum operating pressure, a cryogenic tunnel uses less total energy, and 
therefore costs less to operate, than an ambient temperature tunnel doing the 
same amount of testing. Even so, the operation of large cryogenic tunnels will 
be very expensive in absolute terms. For example, the NTF, when operating at 
its maximum Reynolds number of 120 x lo® at transonic speeds, will use liquid 
nitrogen (LN 2 ) at the rate of approximately 454 kg/sec (1000 Ib/sec) . Although 
the average IW 2 usage rate in the NTF will be much less than 454 kg/sec, it is 
still highly desirable, if not essential, that the tunnel be operated in an 
efficient automatic manner in order to minimize LN 2 consumption. 

Experience with the Langley 0.3-Meter Transonic Cryogenic Tunnel (0.3-m 
TCT) (ref. 2) demonstrated the need for automatic control systems if the full 
potential of the cryogenic tunnel is to be realized. Although the basic cryo- 
genic process is stable, there is a strong tendency for drifting from preset 
test conditions because of the heat-transfer effects created by imperfect tun- 
nel insulation. Additionally, due to the high degree of process cross-coupling 
between independent control variables (temperature, pressure, and drive-fan 
speed) and the desired test conditions (Mach number, Reynolds number, and 
dynamic pressure) , man-in-the-loop operation has proven to be much less effi- 
cient than desired in terms of LN 2 3^*3 electrical-power usage. 

In order to meet the specific control needs of the 0.3-m TCT as well as 
to study the control requirements of cryogenic tunnels in general, a study has 
been undertaken at LaRC to develop and validate a mathematical model of 0.3-m 
TCT process, to utilize the model in a hybrid-computer simulation to design 
temperature- and pressure-feedback control laws, and to evaluate the adequacy 
of these control schemes by analysis of closed-loop data. Hie results of this 
study are reported herein. 


SYMBOLS 

A area, m^ or percent 

b pressure-loss coefficient 

Cj„ specific heat capacity of metal, J/kg-K 

Cp specific heat capacity of gas at constant pressure, J/kg-K 

Cy specific heat capacity of gas at constant volume, J/kg-K 

D dimensionality constant 

E total energy, J 

h specific enthalpy, J/kg 
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Subscripts 

a 

c 

e 


= 597(1 - 0.3M)p-0*035 
Mach number 

mass flow rate, kg/sec 
fan speed, rpm 

total pressure, atm (1 atm = 101.3 kPa) 

heat flow, J/sec 

pressure ratio 

Lapl ace var i abl e 

total temperature, K 

time, sec 

internal energy, J 
specific internal energy, J/kg 
volume, m^ 
mass, kg 

heat-transfer coefficient, J/sec-K 

cooling capacity of gaseous nitrogen, J/kg 

cooling capacity of liquid nitrogen, J/kg 

ratio of specific heats 

fan efficiency 

thermal mass, J/kg 

transport time lag, sec 

acoustic 
tunnel circuit 
exit 
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F fan 

G gas 

i inlet 

L liquid 

m metal 

p plenum 

set set point value 

The use of a dot over a quantity denotes the derivative of the quantity 
with respect to time. 


MATHEMATICAL MODEL OF THE LANGLEY 0.3-METER TRANSONIC 

CRYOGENIC TUNNEL 

The development of a control-ccmpatible model of the 0.3-m TCT (shown 
schematically in fig. 1) has been discussed in detail in reference 3. The 
fundamental principles, assumptions, and equations of this modeling effort 
are presented in the appendix of this paper . These governing equations pro- 
vide a simple, single lumpjed-parameter multivariable process model of the 
0.3-m TCT. Control analysis of the mass-energy interaction existing between 
the mass of gas in the tunnel and the control input parameters was used to 
develop the dynamic model from which the real-time interactive simulator was 
assembled. This resulted in a concept of energy control of the tunnel by 
manipulation of the process states. 

The basic thermodynamic equation for determining the time rate of change 
of the total energy of the gas (ignoring pxstential and kinetic energy) is 


E = Heat fran the tunnel metal + Heat of compression from the fan 

+ Heat frcsn injection of LN 2 + Heat from exhaust of gaseous nitrogen (GN 2 ) 

or 

g p Surface 

E = — (WqCvT) = -J Qm + ^ + niL^^L “ ^°pT (D 

with heat energy added to the gas being positive. 


4 




<■ V' 


After mathematical manipulation, adding the transport delays involved in 

the measurements, and using the definitions Q™ = , ot = (cn - c„)T, and 

1 + tms ^ ^ 

3 = h^ - CpT from reference 3, the temperature dynamics of the tunnel can be 
wr itten 


3t/ 1 + tQs\ 
3tV + t„,s/ 


(B + a) 
6 


m^e 






( 2 ) 


where 6 = WqC^ + WmCm* 

From the ideal gas law, the pressure of a confined gas is proportional to 
the mass of the gas and its temperature, p = K.iWqT, After taking the appro- 
priate derivatives and making substitutions, the time rate of change of pressure 
may be expressed as 


9p p 3 t p 

3t T 3t Wq 3t 


After substitutions and manipulation, the pressure dynamics for the tunnel can 
be written 


P = 


op 


3t Wg 


m^e 


L T O ^ 

^ mce 

Wg 






3t 


(4) 


/ 7pn\ 

where b = 0,19711 - I, Equations (2) and (4) permit the description of the 

temperature and pressure response characteristics of the gas in the tunnel due 
to the tunnel control inputs. 


Other important relationships necessary for modeling the 0.3-m TCT are 
given as follows. The form of the equations used here reflects the specific 
physical characteristics of the 0.3-m TCT and the various elements used for 
measurement and control of the test conditions. 


Mach number : 


M 


Ne 


s 


Km'/T( 1 + tps) 


(5) 
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vrtiere 


Kjj = 597(1 - 0.3M)p-0-035 

tp plenum time constant 

Tg acoustic time lag 

Tunnel circuit time; 


^c “ 


0,0149V/ 


AM\/t \ 


jl + 250 — j(l + 0.2 m2)' 


( 6 ) 


LN 2 flow rate; 

= Kl^L = 3.47\/Pl - p Al (0 ^ Al ^ 100%) (7) 

GN 2 flow rate; (Since Kq is a function of p, two relationships are appli- 
cable for GN 2 flow rate.) 


p p 

me - Kq — Aq = 2,18 — Aq 

VT VT 


p > 1 . 5 atm 
0 ^ Aq ^ 100% 


p 

/I .sy -7 

II 

ll 

II 


VT 

L V H J 


21 , 8 Aq 

{i 


^ p ^ 1 . 5 atm\ 
^ 0 ^ Aq ^ 100% ) 


( 8 ) 


Mass of gas in tunnel; 


Wq = 341.4 —fl + 250 — Vl - 0.033m'' -S) 

T V t2/ 


(9) 


Fan pressure ratio; 


r 


0.197 1 



(10) 
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Metal enthalpy; 


h = 2,75 t2 - 0.0026t3 


( 11 ) 


Fan dynamics; 


N 1 

Nset 0.2 s 2 + 0.56s + 1 


( 12 ) 


Cooling capacities; 


6 = - 1 21 + 2p + 


66 


391 p 


+ 1 . 04T ah Cr,T 


Ot — (^p “ Cy) T 


Fan heat; 


Qp = 


KppM^ \Jt 
(1 + 0.2m2)^ 


(13) 


(14) 


where 


Kp 


A 


6965ACpb 

n 



These equations describe the lumped-parameter mathematical model of the 
0.3-m TCT shown in the sketch in figure 1 and the photograph in figure 2. The 
location of the three control inputs, the sensor locations, and some of the 
physical characteristics of the tunnel are shown in figure 1. The multivari- 
able model used for developing the simulator is presented in figure 3 in matrix 
form. The model outputs (temperature rate, Mach number, and pressure rate) are 
related to the inputs (I 1 N 2 injection-valve area, fan speed, and GN 2 exhaust- 
valve area) through the complex dynamics of the system. Examination of the 
model shows it to be highly coupled and nonlinear in nature. If coupling 
between the variables were not present, that is, if there were no off-diagonal 
matrix terms, variation of the EN 2 injection-valve area, the drive-fan speed. 
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or the GN 2 exhaust-valve area would result in independent changes in the total 
temperature, test-section Mach number, and total pressure, respectively. 

The valve-area openings in the matrix control the mass flow rates of liquid 
and gaseous nitrogen through the process dynamics to affect rates of change in 
temperature and pressure. First-order changes in Mach number are controlled by 
changing the speed of the drive fan. An important feature of this model is the 
large contribution to the process dynamics made by the energy stored in the 
metal of the tunnel and by the GN 2 heat-transfer time constant (which varies 
significantly as a function of test-section Mach number, pressure, and temper- 
ature) . The importance of these factors will become evident during the exam- 
ination of tunnel response data in a subsequent section of this paper. 

Transport time delays have been included in the model to account for the 
mean time between initiating the various control inputs and measuring the 
resultant changes in the test conditions. For example, the transport time delay 
associated with the input of I 1 N 2 is the lag corresponding to the transit time 
from the IiN 2 injection valves to just ahead of the screen section where temper- 
ature is measured by a thermocouple. 


Simulation of the 0.3-m TCT 

At the beginning of this study, it was decided to develop a computer simu- 
lator of the 0.3-m TCT to give hands-on interactive capability for the process 
and control studies. Furthermore, it was anticipated that a simulator of this 
type would be valuable for evaluating various operating techniques as well as 
for the training of tunnel operator personnel. Consequently, the computer 
selected for the modeling and simulation effort was a hybrid computer system 
which permits manual interaction along with high-speed digital computation. 

Solution of the system process equations was undertaken on the hybrid- 
computer system shown in figure 4. Basically, the computer consists of a cen- 
tral digital processor unit and an analog computer electrically interfaced for 
data exchange. The rate of update between the digital and analog computer sys- 
tems is 25 Hz. 

Because of the nonlinearity and the time dependency of the model, compu- 
tations were segregated into two parts. All nonlinear computations were per- 
formed digitally but integration of process equations took place continuously 
in the analog side of the computer system. Shown in figure 5 is a functional 
block diagram of the hybrid-computer simulation. Ten inputs from the analog 
computer unit (after analog-to-digital conversion) are used in the digital 
computer to compute the equations shown in figure 5. The computation in the 
digital computer results in 11 outputs appropriately scaled and converted to 
analog signals for input to the integrating analog side of the hybrid system. 
The loop diagrams of the model dynamics and simulation control along with the 
digitally controlled attenuators which allow the nonlinear digital calculations 
to be introduced into the analog computer are presented in figure 6. Tunnel 
temperature dynamics (fig. 6(a)) are obtained from integration of the temper- 
ature rate whereas tunnel gas temperature T is used to derive the metal 
temperature by using the inverse of the metal time constant l/t^. The metal- 


8 



if 


to-gas heat flow signal Tm - T is used to derive -T in a looped manner. 
Measured gas temperature for display and closed-loop controller use is created 
by digital data provided through the digital attenuators. 

Pressure, fan, and Mach number dynamics (fig. 6(b)) are derived by the 
integration of scaled pressure rate lOp. A coupling term N has been included 
to generate the influence of fan acceleration on the rate of change of pressure. 
Fan-speed modeling is effected by an acceleration-limited second-order transfer 
function. Fan revolutions per minute and Mach number display information for 
the simulator operator panel are acquired after appropriate scaling with analog 
potentiometers. 

Figure 6(c) shows details of the patching necessary to create and control 
the 3-sec pulse inputs used to study the tunnel-model behavior from balanced 
conditions after various disturbances. Also shown in this figure are the valve- 
area display signals for the operator console panel, which will be described 
later . 

The last analog diagram used by the hybrid-computer simulator and shown 
in figure 6(d) is the closed-loop feedback control law for the temperature and 
pressure loops. Proportional and integral (PI) type controllers and the gain 
schedules for the respective variables are displayed. Voltage clipping was 
necessary to prevent negative voltages since the valve-opening area can only 
be positive. 

It was recognized early in this work that experimental data were required 
to validate the process model. As a result, the first phase of work was to 
obtain data for dynamic response of temperature, pressure, and Mach number to 
changes of the input variables throughout the tunnel operational envelope. 

Figure 7 is a photograph of the control panel of the 0.3-m TCT as it was at 
the time of these response tests. Principal control inputs are IiN 2 injected 
rate, GN 2 exhaust rate, and fan speed. Both the LN 2 and GN 2 flow rates are 
manually controlled by commands to specially designed digital valves. These 
digital valves are revolutionary advances in the state of the art of fluid con- 
trol and measurement. Multiple-venting control ports, each fitted with an 
accurately precalibrated, binary-weighted venturi, allow selection of 256 
or more discrete flow rates. Thus, for the 8-bit digital valves used in the 
0.3-m TCT, fluid flow rate can be controlled to within 1 part in 256 (2®) with 
extreme repeatability. Because each control port is bistable (either fully 
open or fully closed) , there is no valve deadband or overshoot which can be 
experienced with conventional analog valves. Also, each port can be modu- 
lar ly replaced or serviced without removing the whole valve from its location. 
Figure 8 is a schematic of the commercially available digital exhaust valve. 

Inputs to the injection or exhaust systems were realized by pulsing a 
selected element of the appropriate digital valve. For these tests, the 
6.25 percent of full-flow element was cycled. Mach number dynamics were 
determined by means of a 100-rpm step change in speed of the drive-fan motor. 
Total temperature, total pressure, and Mach number as well as the three input 
commands were recorded at a rate of 40 samples per second for the tunnel con- 
ditions (given in table I) , which span the entire operational range of the 
0.3-m TCT, The photograph of the tunnel (fig. 2) shows the I 1 N 2 injection 
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valves (upper right), the fan section (lower right), and the piping leading to 
the digital exhaust valve (center left) . 


Operator Control and Display Panel 

In addition to studies of the process dynamic behavior, an important 
reason for the simulation was to achieve man-in-the-loop interaction so that 
0.3-m TCT operational techniques could be investigated. Also, with such a 
simulator operator training could be accomplished at an insignificant cost 
compared with training in the tunnel. After appropriate scaling and electri- 
cal hookup to the analog patchboard, an operator display and control panel 
(shown in fig. 9) was made available for the hybrid-computer simulation. 
Operator inputs are made through potentiometers which control the simulated 
liquid- and gaseous-nitrogen valve openings and the fan speed. The temper- 
ature, pressure, and Mach number digital displays permit realistic simulation 
of tunnel operation including operator imposition of constraints such as maxi- 
mum cool-down rate or metal-to-gas temperature difference. Displays not pro- 
vided on the simulator control panel but present on the tunnel control panel 
include liquid-nitrogen pump pressure, drive-fan motor power, and digital- 
valve status indicator lights. 


Recording Instrumentation 

The technique used to record the experimental data acquired from the 
0.3-m TCT was straightforward and consisted of simply patching analog param- 
eter outputs to a remotely located recording system. Data were directly 
digitized at a rate of 40 samples per second, converted into engineering units, 
and placed on computer files so that computer plots could be generated and 
analyzed relative to the simulator results. 

The procedure for obtaining the simulator data was not quite as straight- 
forward, Instead of direct digitization, the simulation data had to be first 
recorded on an FM wide-band recording tape. The recorder is shown in figure 10. 
The data were then put through a 40-sampler-per-second digitizing process for 
conversion to engineering-unit tapes which were used to create the computer 
file data base for simulator validation. 


Validation of the 0.3-m TCT Tunnel Simulator 

In order to validate the analytical model of the 0.3-m TCT used to develop 
the real-time cryogenic-tunnel simulator, a number of tests were performed both 
on the tunnel and the simulator. The results of these separate tests are com- 
pared and discussed. Details of the experimental-response tests are presented 
and the features of the data are correlated with the mathematical terms for 
validation of the simulator model. 

It is well known that linear time-invariant dynamic systems yield tran- 
sient characteristic responses to deterministic inputs (i.e., sine-wave, step. 
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ramp, impulse, or doublet inputs) . These time-dependent responses can be 
directly compared with simulated responses. The cryogenic process is obviously 
time invariant, and though it is nonlinear, the tunnel characteristics can be 
easily analyzed by assuming linear behavior for small perturbations around an 
equilibrium point. Because the cryogenic- tunnel process has open-loop integrat- 
ing temperature and pressure dynamics, step-input disturbances from equilibrium 
were deemed unsatisfactory since saturation of responses would result. Also, 
for the actual tunnel, sinusoidal disturbances are difficult to realize except 
for fan-speed inputs. Consequently, easily generated mass impulse inputs were 
selected for creation of transient responses for both the 0.3-m TCT and the 
simulator. An additional argument for impulse inputs is that responses can be 
analyzed with regression analysis techniques for "parameter identification even 
when the responses are corrupted by measurement noise. 

During the response tests, both the 0.3-m TCT and the 0.3-m TCT simulator 
were established at specified equilibrium conditions by open-loop adjustment 
of the three control inputs and then subjected to impulse disturbances of LN 2 
mass injection, GN 2 mass exhaust, and drive-fan speed. The disturbances were 
applied serially with sufficient time allowed between impulses to allow tran- 
sients to decay and equilibrium to be reestablished. 

The amplitudes and duration of the perturbations were selected so that 
changes in the tunnel conditions would be in excess of system noise. Thus, 
a 6.25-percent increment in the LN 2 valve full-open area, a 6,25-percent incre- 
ment in the GN 2 valve full-open area, and a 100-rpm decrease in drive-fan speed 
were chosen as the input perturbation amplitudes. Each of these inputs was 
imposed for a period of approximately 3 sec. This pulse duration corresponds 
to 3/4 to 4 tunnel-circuit time periods, depending on test conditions. Typi- 
cally, these pulse inputs create changes in temperature of 0,5 K to 6.0 K, 
changes in pressure of 0,017 to 0.102 atm, and changes in Mach number of 0,006 
to 0,030. The dynamic responses of the experimental data and the simulation 
data are compared in figures 11 to 37 for the various test conditions selected 
to cover the full range of operation of the 0,3-m TCT. 

Inspection of the data shows that the records for the simulator responses 
are practically noise free whereas the tunnel data exhibit some measurement and 
recording noise. The Mach number traces show a larger level of noise because 
Mach number is calculated from two pressure records, each of which has its own 
noise component. One particular problem encountered in obtaining the experi- 
mental data from the 0.3-m TCT was a 10- to 20-percent random drift from the 
nominal value in the IiN 2 pump pressure. This prevented accurate measurement 
of the mass flow of the LN 2 injected into the tunnel. 

Once equilibrium was established, perturbation of the 6.25-percent element 
of the IW 2 valve admitted an additional mass of LN 2 into the tunnel circuit, 
typically between 0.5 and 1,5 kg per pulse. Perturbation of the 6.25-percent 
element of the GN 2 valve exhausted an additional mass of GN 2 to the atmosphere, 
typically between 0.2 and 1 , 0 kg per pulse. 
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ANALYSIS OF TUNNEL AND SIMULATOR 


Transient Responses 

Figures 11 , 28 , and 29 are provided as typical impulse responses alpha- 
betically subdivided so that various features and characteristics of the 
temperature, pressure, and Mach number trajectories can be identified and 
explained. Prior to any intentional system disturbance the equilibrium con- 
ditions for the cases studied were established. The segment ABCD defines the 
portion of the response influenced by the perturbation to the LN 2 valve area. 
The segment EFG defines the portion of the response influenced by the pertur- 
bation to the GN 2 valve area. The segment HJK defines the portion of the 
response influenced by the perturbation of the drive-fan speed. Correlation 
of these dynamic records with the model mathematics was performed to provide 
better comprehension and insight into why specific results are obtained. This 
is done in the following sections for the three control inputs. 


LN 2 Pulse Input 

Temperature response .- Refer to the temperature trace of figure 11. Prior 
to A the tunnel was placed in a balanced mass and energy condition such that 
the rates of change of pressure and temperature were zero (p = T = 0) , the 
liquid-mass inflow was equal to the gaseous-exhaust outflow ( 111 ^ = m^) , and 
temperature, pressure, and Mach number were constant. During segment AC, LN 2 
mass flow rate is increased by Am^ due to the increase in valve area of 
6.25 percent of the full-open valve area during the perturbation. This creates 
a negative temperature gradient (T < 0) due to the negative 3 term in equa- 
tion (2) . As the temperature of the gas in the tunnel rapidly drops, the tem- 
perature of the metal in the tunnel shell lags behind and begins releasing 
stored heat energy to the gas. Eventually the amplitude of the temperature 
gradient is dynamically reduced as dictated by the lead/lag term of equa- 
tion (2) . Gas temperature (at the measuring station) continues to drop beyond 
point C. The value Am^ returns to zero due to removal of the perturbation 
and closure of the 6.25-percent element of the LN 2 injection valve. Recovery 
from the perturbation takes place in the segment CD due to continued heat trans- 
fer from the walls to the gas. At point D, most of this heat transfer is com- 
pleted. It may be noted that there exists an obvious temperature transport 
delay relative to the commanded input position. This corresponds to that por- 
tion of the circuit time constant (see eq. (6)) consistent with the circuit 
length between the LN 2 injection station and the temperature-measurement sta- 
tion. All temperature data (figs. 11 to 37) display this delay characteristic. 

The temperature differential established between C and D by first a neg- 
ative gradient (due to liquid input) and then a positive gradient (due to wall 
heat release) can be predicted by reexamination of equation (2) . For zero 
change in exhaust flow and no added heat contribution from a steadily running 
fan, only LN 2 influences the differential temperature AT (or T) . If AT 

(1 + t„,s)/3 +a\. _T-_s l/l + tn,s\ 

is represented by mre ^ and - is transiently 

(1 + tQs) V 0 / 0 \1 + tQsJ 
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approximated by 


1 


(since the energy of the 


term in 


0 


WQCy 


is released only after a long time period). 




then 


(BmLe + Ql) 


Wg^v + WmCjn 


By using 


WqCy 

this simplified analysis, the largest value of At should occur when 3 is 
largest and Wq is smallest. Equation (9) shows Wq is smallest when pres- 
sure is lowest and temperature is highest. For the tunnel conditions tested, 
Wq is smallest at a pressure of 1.5 atm and a temperature of 275 K. The 
experimental (tunnel) data presented in figures 11 and 29, for example, con- 
firm that the maximum AT does, in fact, occur under conditions of minimum 
Wq. For the condition of constant temperature, as the mass of the gas in the 
tunnel increases with increasing pressure, the amplitude of At should 
decrease correspondingly. Such a decrease in At is shown, for example, in 
figures 29, 32, and 35. Thus, the smallest At occurs when 3 is smallest 
and Wq is largest. For the conditions examined this occurs at 5.0 atm and 
100 K. The tunnel data presented in figures 17 to 19 exhibited the smallest 
values of At as predicted. The magnitude of recorded temperature differ- 
entials is somewhat smaller than predicted because of some filtering of the 
actual temperature measurements due to the response character istics of the 
thermocouple used for the temperature sensor . 


Pressure response .- Refer again to figure 11. With attention now focused 
on the pressure trace, the following analysis will explain the pressure behavior 
of the tunnel. Equation (4) is the tunnel state relationship for the rate of 
change of operating pressure. With constant values of gaseous exhaust and fan 


contributions in the pressure rate equation, only the cooling term - T and 

T 

P • — Tt s 

the mass term — Amj-e ^ influence the pressure gradient. For segment AB 

Wg 


of the pressure trace, p is negative because the cooling term is larger than 
the mass term. At point B the two terms are equal and opposite. As a result, 
p = 0 and the minimum pressure is realized. For the segment BC, the mass term 
becomes much_ larger because of the additional liquid-mass input, and despite 
a negative T, a positive p ensues. Thus, the positive pressure gradient 
from B to C is due to the dominant heat transfer from the tunnel wall to the 
test gas. At point C, Ahil is reduced to zero by reducing the valve incre- 
mental area to zero. Because the wall is warmer than the gas, the net heat 


p . 

flow to the gas causes T to become positive in the - T term. As a result, 

T 

p increases until the heat transfer is complete and a new equilibrium value 
of p is established. 


This initial decrease in pressure due to the liquid input followed by an 
increase in pressure is characteristic of all conditions tested for the simu- 
lator and the tunnel. (See figs. 11 to 37.) The location of the minimum pres- 
sure (point B) varies as a function of the test conditions of the tunnel. For 
gas temperatures of 100 K, the minimum pressure consistently occurs very near 
the end of the liquid input (point C) . (See figs. 11 to 19.) On the other 
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hand, at a gas temperature of 275 K, the minimum pressure is very near the 
beginning of the liquid input (point A) . This indicates that a much more rapid 
pressure recovery is established because of the large amount of energy released 
from the tunnel walls. This is as expected since the metal enthalpy (eq. (11)) 
at 275 K is eight times the enthalpy at 100 K. At the intermediate gas temper- 
ature of 200 K, the minimum pressure varies in location between points A and C. 
This is due to the different values of the cooling effect experienced because 
of varying degrees of wall heat release. It should be noted that minimum pres- 
sure shifts toward point C as gas-weight increases (i.e., increase in pressure) . 
The data presented in figures 20, 23, and 26 illustrate this for the 200 K case. 

An interesting point in the study of the response of tunnel pressure to 
liquid input is the system behavior at point D. When this point is reached, 
most of the heat-transfer effects have been completed and tunnel recovery 
begins. At temperatures of 200 K or 275 K (figs. 20 to 37), the pressure level 
at point D is much higher than the starting or initial equilibrium pressure 
whereas at 100 K (figs. 11 to 19) the final pressure is almost the same as the 
starting pressure. On the other hand, the recovery pressure gradient is not as 
pronounced at the lower temperature. The higher temperature data indicate that 
the slope of the recovery gradient increases more negatively as the Mach number 
increases. (See figs. 20 to 22 and figs. 29 to 31.) These larger negative 
gradients can be explained with the help of equations (4) and (8) . The rate of 
change of pressure p is a function of exhaust-mass flow rate which itself is 

P 

related to — Aq. The exhaust valve area Aq at the higher Mach numbers has 

'It 

to be larger to maintain equilibrium, and therefore when the incremental liquid 
input is removed, more GN 2 is being exhausted from the tunnel, resulting in 
larger negative values of p. Note that there exists essentially no delay in 
the response of tunnel pressure to changes in IN 2 flow rate. 

Mach number response .- The response of Mach number to liquid input can be 
evaluated by using equation (5) , which indicates that Mach number is inversely 
proportional to the square root of temperature for constant fan speed N. 
Therefore, when the IN 2 flow rate is increased, temperature falls, the speed of 
sound decreases, and since flow velocity is constant, the Mach number increases. 
This situation is illustrated, for example, in figure 11 in the increase in Mach 
number for segment AC. Coupling between temperature and Mach number is caused 
by the increase in fan pressure ratio due to the reduced temperature. (See 
eq. (10).) At point C, the Mach number response has the largest deviation. 

Am At 

which can be related to the gas temperature-response amplitude by — = - — 

M 2T 

(derived from eq. (5)). Figure 11 shows that for a Mach number of 0.3 and a 
temperature of 1 00 K the perturbation in liquid flow rate causes a change in 
Mach number of 0.006. (A computed value of At = -4 K results.) Because of 
the previously mentioned first-order filtering effects of the thermocouple, the 
recorded experimental data from the tunnel indicate a At of only about -1.5 K. 
All of the temperature-response data from the tunnel presented in figures 11 
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to 37 are subject to this type of filtering. However, there is generally good 
agreement between tunnel and simulator responses. 

In summary, test-section Mach number and pressure responses of the tun- 
nel to an LN 2 input agree quite well with the mathematical model predictions 
described by the simulator responses for the entire operational envelope. 
Actual tunnel temperature responses did not duplicate simulator predictions 
quite so well due to the dynamic characteristics of the thermocouple used for 
temperature measurements. 


GN 2 Exhaust-Pulse Inpu.t 


Once again in reference to figure 11, the tunnel and simulator responses 
are examined for an exhaust pulse of about 3 sec. These responses are marked 
EFG for selected study. Equilibrium conditions have been established prior to 
point E(p = T = 0, niL = Ag) . As with the perturbation to the LN 2 flow rate, 
the perturbation to the exhaust-gas flow rate is induced by increasing the valve 
area by an incremental 6.25 percent of the full-open value. This was done 
physically for the tunnel and in simulation on the computer for all conditions 
in table I. In the perturbation process an increase in the exhaust-gas mass 
flow AniQ is produced. As described in equation (4) , this results in a neg- 
ative pressure gradient (p = - — AmG® Because the exhausting causes a 

\ Wg j 

negative temperature gradient (T) during segment EF, there is a small second- 


order negative pressure-gradient contribution due to the 


P • 

- T 
T 


term in equa- 


tion (4) , By and large however, the pulse causes a linear pressure drop during 
segment EF. The linear pressure drop consistently repeats for both the tunnel 
and simulator responses as shown in figures 11 to 37. At point F the incre- 
mental valve area is made zero, thereby making Aitiq = 0. Once Aitiq is made 
zero, the slope of the pressure gradient p is characterized by the steady- 
state gaseous-mass flow rate mg, which is a function of p/^. Thus the largest 
pressure gradients exist when pressure and temperature are maximum (p = 5.0 atm 
and T = 275 K) . Examination of the data shows this to be true. In contrast, 
the smallest pressure gradients occur when p = 1,5 atm and T = 100 K. During 
the interval FG, the pressure gradient is a direct function of the gaseous-valve- 
area opening Ag- During steady state, this valve-area opening is largest at 
high Mach numbers in order to maintain equilibrium operation. Because of the 
large opening at high Mach numbers, more GN 2 is being exhausted. When the 
6.25-percent incremental valve-area perturbation is removed, the reduction in 
AniQ is greatest at high Mach numbers. For this reason, the pressure recovery 
during segment FG is fastest at M = 0.75 and M = 0.90 and slowest at 
M = 0.30. 


As expected, in general the perturbations of the exhaust-mass flow rate 
have more effect on pressure than on either temperature or Mach number. During 
segment EF, however, a negative T, which is filtered by the lead/lag term, is 
found. Equation (2) is the governing equation for this temperature gradient. 
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when AniQ is made zero at point F, T becomes positive due to the reduced gas 
mass exhaustion, thereby establishing a new temperature. The magnitude of At 

^ * “T#^S 

can be transiently approximated by mge since there is zero change in 

WqCv 

liquid flow rate and no added heat contribution from a steadily running fan 
(eq. (2)). since iiiq is a function of p/'jT (eq. (8)), AT is expected to 
be largest at high pressures and low temperatures. Comparison of figure 11 
with figure 29 illustrates this effect quite well. 


Fan-Speed Pulse Input 

The trajectories marked HJK in figures 11, 28, and 29 correspond to the 
response of total-pressure, total-temperature, and test-section Mach number 
dynamics of the cryogenic tunnel to a pertubation in fan speed. As with the 
other types of pertubations, system equilibrium is established and during seg- 
ment HJ the fan speed is reduced by 100 rpm. The fan speed decelerates almost 
linearly as described by the fan dynamics of equation (12). Reduced fan speed 
results in a pressure-ratio reduction which manifests itself as a decrease in 
Mach number. Because the fan is accelerated during segment JK back to the 
original speed, the pressure ratio is increased causing Mach number to return 
to its original value. The maximum Mach number deviation occurs at point J. 

An 

This is related to fan speed by Am = — as indicated in equation (5) , Since 

An is made constant for all runs, the largest Mach number deviation takes 
place at the lower temperatures. Comparison of figures 11 to 1 9 with fig- 
ures 29 to 37 confirms this phenomenon. At 100 K, a typical Mach number change 
is about 0.034 compared with a 0.020 change at 275 K. Reduced fan heat produc- 
tion due to the speed decrease contributes a secondary influence on the Mach 
number reduction. The increase in Mach number caused by decreased heat of com- 
pression is quite small and cannot be easily separated from the primary effect 
of the decrease in fan speed. Total- temperature data reflect this by a slight 
temperature decrease during segment HJ and an increase during segment JK. 
Equation (2) indicates the influence of fan heat Qp on the temperature gra- 

KppM^# 

dient. Because fan heat (eq. (14)) is a contributor to equa- 

(1 + 0.2m2)3 

tion (2) , the maximum At occurs when p, T, and M are maximum. A com- 
parison of figure 37 (5,0 atm, 275 K, and M = 0.75) with figure 11 (1.5 atm, 
100 K, and M = 0.3) demonstrates this clearly. The transit-time delay previ- 
ously mentioned is also noticeable in the temperature traces of these records. 

Except for figures 12, 25, and 29, the tunnel and simulator responses 
agree quite well. In figure 12, the oscillatory pressure and temperature 
responses for the tunnel are due to an inadvertent cycling of one of the 
smaller exhaust-gas valve elements which is not shown in the figure. Two 
fan-input commands were made which caused cumulative effects in the tunnel 
variables shown in figure 25, The second input, which was inadvertently made, 
is not shown in the figure in order to show only the effects of the inten- 
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tional input. The large pressure and temperature magnitudes experienced in 
figure 29 are a result of a 5-sec liquid pulse rather than a 3-sec input. 


Quasi-Steady-State Warm-Up and Cool-Down Tests 

In addition to the transient-response studies/ tests on cool-down and 
warm-up of the tunnel with excess IiN 2 flow and excess fan heating, respectively, 
were conducted to validate the quasi-steady-state heat-transfer model. 

The basic energy equation of the cryogenic-tunnel process is 


dT 

WqC« t Cl 
dt 


^Surface 

^T(mL - me) = - \ Qm + Qf + " "kSCpT 


(15) 


where 


/^Surface Wj^CjnTs 

J Qn, = r— = (T - T^)y 


1 + tmS 


In order to validate this heat-transfer model, time histories of the average 
metal temperature Tm and gas temperature T were recorded for ordered input 
energy terms. A match between temperature, temperature difference, and time 
was sought for the tunnel and simulator quasi-steady-state responses. 

Cool-down .- The actual cool-down process used for the 0.3-m TCT consists 
of injecting liquid nitrogen at a constant rate and in a sufficient quantity 
so that the cooling capacity of the liquid is greater than the heat conducted 
through the tunnel shell and the heat generated by the fan [|mL(hL - CpT) » Qp]| . 
This condition is easily realized by running the tunnel at a very low Mach 
number while injecting the LN 2 . If we neglect the heat conducted through the 
tunnel walls, this implies that the energy equation can be simplified to 


pSurface 

WqCv — “ " \ ^ (1®) 

dt 


Both the tunnel and the simulator were brought to an equilibrium condition of 
290 K and 1.5 atm while maintaining a constant fan speed of 1200 rpm. Liquid- 
nitrogen injection pressure was held constant at 6 atm throughout the cool- 
down. Starting from a condition of equilibrium, the actual cool-down process 
was started by increasing the LN 2 injection-valve area to 12.5 percent to 
achieve a 1.02-kg/sec liquid-nitrogen flow rate. Tunnel pressure was main- 
tained constant during the cool-down process by manipulation of the exhaust- 
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valve area. The temperature of the gas, the average temperature of the metal 
of the tunnel shell, and the test-section Mach number were recorded during the 
cool-down. Both metal and gas temperatures were reduced and cool-down was con- 
sidered complete when the gas reached a tCTiperature of about 105 K. The total 
amount of LN 2 consumed was also determined and recorded. Figure 38 shows the 
cool-down profile for both the tunnel and the simulator. As can be seen, both 
gas and metal temperature trajectories match quite well. An initial gas-to- 
metal temperature difference of 65 K exists. This difference decreases during 
the cool-down, with a difference at the end of the cool-down process of only 
15 K. The I 1 N 2 consumption to cool the tunnel from 290 K to 100 K was about 
1825 kg, whereas the computer simulation for the same cool-down predicted an LN 2 
usage of 1800 kg. This close agreement between the simulator prediction and 
the tunnel data gives confidence that the quasi-steady-state heat- transfer 
model is accurate. 

Warm- up . - For a typical tunnel warm-up, the LN 2 flow is reduced to zero 
while the tunnel is operated at constant Mach number and pressure. For this 
study, the warm-up was made after establishing equilibrium conditions with both 
metal and gas temperatures at 100 K. Since the LN 2 flow rate is zero and GN 2 
flow rate is very nearly so, the heat generated by the drive fan is the dominant 
factor in the warm-up process. Neglecting any heat addition through the tunnel 
walls, this results in the simplified energy equation 


(3T ^Surface 

WqCv - ^ ^ J 2m H7) 


where 







(T - Tm)y 


and 


KppM^ \/t 

Of = —2 >> "teCpT 

(1 + 0.2m2)-^ 


Both tunnel and simulator were set at equilibrium conditions of 1 00 K and 
2 atm. Mach numbers of 0.5, 0.6, and 0.7 were used for the computer simulation 
in order to show the effect of Mach number on the warm-up profile. Once equi- 
librium was established, the LN 2 flow was reduced to zero by closing the LN 2 
injection valve. Mach number and pressure were maintained constant throughout 
each warm-up. Figure 39 shows the warm-up profiles predicted by the simulator 
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for the three different Mach numbers. As expected, the data indicate larger 
gas-to-metal temperature differences as well as faster warm-up for higher Mach 
numbers. Figure 40 shows a comparison of an actual 0.3-m TCT warm-up profile 
and a corresponding simulator warm-up profile at 2 atm and M = 0.6. As can 
be seen, the tunnel warms up slightly slower than predicted by the simulator, 
but the metal-to-gas temperature differences between the tunnel and simulator 
agree to within 2 K. These differences between the warm-up profiles are attrib- 
uted to slight variations in tunnel pressure at the beginning of the warm-up 
process and perhaps a small error in the constant Kp. 

Generally, the behavior of the 0.3-m TCT and the behavior predicted by the 
computer simulation of the tunnel show good agreement not only in the transient 
responses, but also in the quasi-steady-state cool-down and warm-up character- 
istics. This agreement has been found over the entire operational envelope. 

In view of this agreement, the mathematical model developed to describe the 
tunnel is assumed to be valid and acceptable for transient, quasi-steady-state, 
and steady-state performance. 


Model Closed-Loop Analysis 

Since the beginning of operation of the 0.3-m TCT, the need for improving 
the quality of data obtained during tunnel tests has become increasingly appar- 
ent. A major improvement in the quality of the data can be realized if the 
test parameters can be held to closer tolerances for longer periods of time. 

It is reasoned that if inaccurate and inefficient manual operator manipulation 
for holding and changing the test parameters can be eliminated by the use of 
advanced control schemes, data quality will significantly improve, with a 
simultaneous reduction in operating costs. As shown in reference 4, by reduc- 
ing the run time, one could proportionately reduce the cost of testing because 
of decreased electrical-power and LN 2 consumption. These arguments all alluded 
to the implementation of closed-loop control systems, which are inherently 
faster and more accurate (and generally more reliable) than the human operator. 
Consequently, proportional-integral-derivative (PID) pressure and temperature 
control laws were developed and designed for analytical evaluation of the 
hybrid-computer simulator of the 0.3-m TCT. Figures 41 and 42 show schematic 
diagrams of the pressure and temperature loops, respectively. Shown in both 
figures are the equations describing the digital valves and the respective lags 
associated with the temperature and pressure sensors. Contributions of the 
various input/output elements are shown with appropriate signs conveying the 
influence they have on the process loop. 


The closed- loop mathematical models of the 0.3-m TCT were operationally 
analyzed assuming local linearity and no coupling between controls. For pur- 
poses of illustration, the pressure loop consists of first-order sensor 
measurement dynamics with time constant t 2 » which determines feedback pres- 
sure that is compared with the pressure set-point value. The generated error 
signal drives the PID controller which, after appropriate gain selection, con- 
trols the operation of the nitrogen-gas exhaust value. Root locus analysis of 


the pressure loop indicates that a gain schedule of 
performance . 



provides acceptable 
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A feedback lcx3p similar to that used for pressure control is used for the 
temperature control with the applicable measurement sensor time constant and 
digital-valve dynamics. Included in the temperature control loop is a fan 
power feed-forward term necessary to account for the effect of fan-generated 
heat on the temperature loop. Root locus analysis of the temperature control 

M\|p 

loop indicates that a gain scheduling of — ^ would be satisfactory. 

Figure 43 shows data from the hybrid-computer simulator for the closed- 
loop performance for both low and high tunnel temperatures using only pro- 
portional and integral signals within the pressure and temperature controls. 
Closed-loop data from the simulator were acquired for temperature and pres- 
sure set-point changes and finally for fan-speed changes. 

In figure 43(a), the temperature set point is reduced by approximately 8 K, 
The temperature reaches 82 K in about 10 sec and then slowly drifts downward 
over the next 40 sec. This lag .is due to the heat transfer from the tunnel 
wall. The pressure control maintains pressure to within 0.03 atm and suppresses 
the fluctuations in pressure due to the temperature disturbance in about 20 sec. 
It should be noted that Mach number increases because of the reduction in tem- 
perature. Following stabilization of the gas temperature at 82 K, the temper- 
ature set point is returned to 90 K. The gas temperature reaches 90 K in 5 sec 
and, since the time spent at the lower temperature was too short to reduce the 
metal temperature much below 90 K, the gas temperature stabilizes quickly at 
90 K. Again the pressure settling time is about 20 sec. 

The next disturbance to be evaluated is a pressure set-point change from 
2.13 atm to 1.76 atm. An overshoot in the simulated tunnel-pressure process 
exists but stabilizes quite nicely. Closed-loop temperature oscillations caused 
by the pressure input settle in about 22 sec whereas Mach number follows the 
temperature oscillations in opposite fashion. After the pressure stabilizes 
at 1.76 atm, the pressure set point is returned to 2,13 atm. The relatively 
slow pressure buildup occurs because it is necessary to increase the mass of 
gas in the tunnel. The rate of increase is entirely dependent upon the mass 
flow rate of LN 2 into the tunnel. Closed-loop temperature control for the 
pressure disturbance is maintained to within 0.5 K and quickly returns to the 
temperature set-point value. 

The last input for closed-loop study consists of a change in drive-fan 
speed to force a change in Mach number from 0.75 to 0.65 quickly followed by 
a change in speed to force a change in Mach number from 0.65 to 0.85. Even 
with the large change in Mach number (0.65 to 0.85), very small temperature and 
pressure excursions occur, indicating acceptable performance of the temperature 
and pressure control systems. 

Figure 43(b) illustrates the control-system performance for set-point 
changes at the high temperature of 263 K. Again the pressure and temperature 
loops perform well, with average settling times of about 20 sec. 
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Based on these simulator-generated results, the control laws derived on 
the basis of single- input, single-output analysis provide good closed- loop con- 
trol of pressure and temperature. Coupling effects between process response 
and control input are evident, but converge very well using simple loop closures. 


CONCLUDING REMARKS 

A hybrid-computer simulation of the Langley 0.3-Meter Transonic Cryogenic 
Tunnel has been developed and fully verified. Comparison of simulation and 
experimental transient response data exhibited very good agreement throughout 
the cryogenic-tunnel operational range. Additionally, quasi-steady- state cool- 
down and warm-up profiles have contributed significantly to the model valida- 
tion by virtue of the small differences seen between actual and simulated 
operation. Therefore, it is felt that the single lumped-parameter nonlinear 
multivariable model has been proven to be globally accurate and reliable in 
duplicating tunnel temperature, pressure, and Mach number process dynamics, 

A plethoric compilation of data spanning the tunnel testing range is provided 
as a data base for future reference. 

Proportional- integral control laws designed using the validated mathe- 
matical model have been exercised using the hybrid-computer simulator and have 
been found to be quite satisfactory for control of both the temperature and 
pressure process responses. As a result, these control laws will be imple- 
mented in software for digital microprocessor controllers for temperature and 
pressure regulation of the Langley 0.3-Meter Transonic Cryogenic Tunnel. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
June 23, 1980 


APPENDIX 


THERMODYNAMIC MODEL OF THE LANGLEY 0,3-m TCT 

In this appendix, a dynamic lumped-parameter model is developed for the 
0.3-m TCT circuit. The need for a relatively simple model which can be numer- 
ically simulated and rapidly executed to permit analytical design and inter- 
active usage excludes an analysis based on nonlinear partial differential equa- 
tions of fluid flow. Lumped-parameter techniques are applied to produce a low- 
order model for the dynamics of the process variables at the test-section 
location in the form of differential-delay /ordinary differential equations. 


Temperature Equation 

The first law of thermodynamics applied to this analysis can be written as 


Q + w = E 


This is an energy equation which states that the change in total system energy 
E is equal to the heat Q added to and the work w done by the selected 
system. The energy term consists of the sum of potential energy, kinetic 
energy, and internal energy. Because the change in elevation as the gas circu- 
lates around the tunnel is small, potential-energy contributions were neglected. 
For the tunnel testing conditions operationally used, the kinetic energy of the 
circulating gas was found to be a small percentage of the total system energy 
and was therefore ignored. 

Since no work is performed by the circulating cryogenic gas, this term of 
the general thermodynamic equation is zero. Therefore, it can be seen that heat 
energy Q is the sole source for change in the level of internal energy of the 
gas. Mathematically this can be expressed in special form as 


Q 


3U 

at 


or 


Q 


a (Wqu) 
at 


(Al) 
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where 


u 

internal 

energy, J 

u 

specific 

internal energy, J/kg 

Wg 

mass of 

the circulating nitrogen 


test gas. 


kg 


The left side of the above energy equation is composed of a number of individual 
sources which collectively contribute to the internal energy of the system. For 
example, two-way heat transfer exists between the ambient environment and the 
gas through the imperfectly insulated tunnel metal structure. Heat of compres- 
sion from the tunnel fan is obtained as the test medium is circulated through 
the tunnel. Also, energy is added by injection of LN 2 and exhaust of GN 2 during 
regulation of temperature and pressure conditions within the tunnel. The sign 
convention selected is such that any energy added to the GN 2 tunnel medium was 
considered positive. 


Expanding equation (AT ) by taking the first partial deviative of the right- 
hand side yields 


. 3u 

Q = Wfi — + 

^ 3t 


u 3 Wq 

3t 


(A2) 


Because 3u A Cy 3 t, equation (A2) becomes 


Q 


WqCv 


3t 

— + uWr; 
3t ^ 


or 


Q = WqC^T = uWq 


(A3) 


The term Wq is really the time rate of change of the tunnel resident gas and 
can be expressed as the difference between the injected- and exhausted-nitrogen 
masses, that is. 


Wq = 


(A4) 


Substituting this along with the appropriate metal, fan, IiN 2 , and GN 2 heat 
energy terms for the left side of equation (A3) gives 
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“Qm + ^ + n>L^L " = WIgc^T + u(niL - me) (A5) 


After transposing and collecting terms, this equation is 


mL(hL - u) - me (ho - u) - 4n + ^ = 


(A6) 


However, for a perfect gas (which can be assumed for nitrogen), u = c^T and 
h (3 = CpT. Therefore, the relationship becomes 


m|£^(hL ~ CyT) mgCCpT CyT) Qji^ + Qp — WqC^T 


(A7) 


Expressions which are functions of local pressure, temperature, and Mach_ number 
are now needed for Qm and Qp in equation (A7) . The development of Qp is 
accomplished in another section of this appendix. A detailed analysis of the 
derivation of the metal-to-gas heat transfer 2m based on a numerical estima- 
tion of the heat-transfer coefficient by Bartz is presented in reference 3, 
Results of that work indicate that 




Surface 



7 + t„,s 


(A8) 


By using equation (A8) together with equation (A7) and rearranging we get 


, / ^m*^ \ . 

T WqCv + = niL(hL - c^T) - mcCCp - c^) T + Qp (A9) 

V ■> + W 


By defining 


a = (Cp - Cv)T 
9 = Wqc^ + WijjCm 


and 


f 66 p 

3 = hr - CnT = -121 + 2p + 391 - + 1.04T 

^ \ T T 
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equation (A9) can be expressed as 




' I — itiL (ot + B) — nigOt + 

+ tm® / 


Qf 


or 


/ WQCvtmsX 


^ ^® 

V 0 


= iriL (Qt + 3 ) - ingOt + Qp 


/ 


(AlO) 


This becomes 


. n + tG ®^ 

V ^m® 


a + B\ . /a' 


0 n ■ [eh ^ \eh 


(An ) 


where 


tr. ^ 


WQC^tni 

*0 


and 


Wg = 


pv/ 

341 .4 —h +250 
T \ 



- 0.033m'' • 




Equation (All) is the fundamental temperature relationship for the model with 
the exception of transport delays needed for realism. Because transit time is 
experienced between control input and measurement, pure transport delays are 
incorporated. Liquid nitrogen is injected downstream of the test section in 
the 0.3-m TCT, and since this effect is measured in the settling chamber a time 

delay of e is used. The fan heat of compression occurs after turn num- 

ber two in the tunnel, and since it also is measured in the settling chamber a 

delay of e is needed. A smaller delay, e is used for the gas 
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effects since gas venting is very near the pressure measurement sensor in the 
settling chamber. By applying these delays and performing a bit of algebra, 
equation (AT 1 ) is now 


T 



+ tms\ 
V + tcs/ 




-^lS 



+ 



+ tjnS\ 


+ tQS 


Qpe 


-ToS 


(AT 2) 


The only remaining unknown in equation (A12) is the heat energy generated 
by the fan while circulating the test gas around the tunnel. This contribution 
can be represented by 

Qp = Cpiiip^i At (A13) 

If we assume that the mass flow rate at the test section m* is the same as the 
mass flow rate at the fan inlet (i.e., m* =,n»F,i)» then the fan-inlet mass flow 
rate can be expressed as 


itiF,i = 6965 ^ AM(1 + (AT 4) 


where p is pressure in atmospheres, A is test-section area in square meters, 
and T is temperature in kelvins. This expression can be found in any thermo- 
dynamic text for the calculation of mass flow rate of an ideal gas through an 
isentropic nozzle. 

The temperature change in equation (AT 3) at the fan section is found by 
using the fan pressure-ratio equation for steady, isentropic compressible flow. 
Mathematically this is 


y 

Pe / 

r = — = — (ATS) 

Pi \Ti/ 
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where Y = — r the ratio of specific heats of the gas under consideration. By 
rearranging this, equation (A15) becomes 






Working only with the temperature term in the above we can write 


Ti + AtNY-I 


Because At A Tg - T^, or T^ = Tg - At, equation (AT 7) can be expressed as 


Tg - At + At\y-i 
T p - At / 


which, after adding the At terms becomes 


' Tg 

Tg - At/ 


27 


APPENDIX 


Solving equation (AT 8) for we get 


At = Tg 


which, after assuming a fan efficiency of ti, becomes 


Te 

At = — 
n 


or 


Te 

At = — 
n 


As with any closed-circuit tunnel, the steady-state fan pressure ratio r 
can be expressed as a function of the integrated pressure-low coefficient and 
the normalized test-section Mach number. Mathematically, 

r = 1 + bM^ 


Introducing this into equation (A19), we get 


Te 

At = — 
n 
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Y-1 


The inner term (1 + bM^) ^ can be expanded in a first-order Taylor series 

( T - A 2 

yielding IbM'^. Substitution into equation (A20) gives 



(A21 ) 


Placing equations (A14) and (A21 ) into equation (A13) results in 


Qp = Cj 


_ — 
P -^3 

6965 — AM(1 + 0.2m2) 

■ '')bM2 

_ /t 

J1 V Y / 


Because Tg = T (the total gas temperature), if we let 
and cancel like terms, the above can be written as 


Kp k 


6965 

n 



Qp = KppM-^d + 


0.2m2) ^\[t 


(A22) 


From a control standpoint, equation (A22) should be represented in terms 
of a system input. From a curve fit of tunnel experimental data, a relation- 
ship between the system input of revolutions per minute and Mach number is 
needed. This function was determined to be 


N 

M = (A23) 

KmVt 

where Kjj = 597 (1 - 0. 3M)p~®* Equation (A23) combined with equation (A22) 

yields 
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Qp = — PNm2(1 + 0.2m2)~-^ (A24) 

Km 


Finally, the T system equation of the matrix in figure 3 can be obtained 
by substituting equation (A24) and the following valve mass flow rate equations 
into equation (AT 2): 




where ^ 3 . 47 - P and A^ = Injection valve opening in percent, and 


P 

^ ^ ^ 



/I .5\1 *7 

where Kq = 

2 - 

L \ P / J 


21 .8 


The equations for m^ and mg 


and A (3 = Exhaust valve opening in percent, 
are equations (7) and (8), respectively. 


Mach Number Equation 

As previously mentioned, the fan speed and the test-section Mach number 
are related by the empirical expression 

N = 597M\/t( 1 - 0.3M)p"°*°^^ (A25) 


For any given Mach number, this relationship resulted in a good fit of the 
0.3-m TCT data (within ±10 rpm) , Because of the plenum volume surrounding the 
test section, the Mach number is dynamically influenced. Previous testing indi- 
cated that this influence is of a first-order type effect because of the net 
flow into or out of the plenum. Also, there is a slight delay associated with 
a fan rpm change that is to be interpreted as a test-section Mach number varia- 
tion. From these two facts, equation (A25) becomes 


M 



s 


Km\/t (1 + tpS) 


(A26) 
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where 

% = 597(1 - 0.3M)p-0*035 

tp plenum time constant 

Tg acoustic time lag 

Equation (A26), which is presented in the text as equation (5), is the second 
expression of the model matrix in figure 3 and relates the rpm input to the Mach 
number output. 


Pressure Equations 

The last dynamic model function used to describe the characteristics of 
the 0.3-m TCT is the pressure equation. From the ideal gas law (which nitrogen 
very closely approximates) , pressure is related to temperature according to 


p = KiWqT 


where Wq is the mass of the gas in kilograms. Simply taking the first partial 
derivative of this with respect to time yields 


3p 

_ - KiT — * K,Wo - + M.E. (A27) 


where M.E. denotes manentum effects. But Wq 
equation (A27) to be 


P 

and T 

K-jT 


P 

KiWq 


will cause 


3p p p 3 t 

— + + M.E. 

3t Wq 3t T 3t 


Since 


9Wq . 

= m^ - mg, equation 


(A28) becomes 


3t 


P . 

Wq 


p 3 t 

+ + M.E. 

T 3t 


(A28) 
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or 


3pp. p . p . 

= — niT — — n\Q + — T + M.E. (A29) 

3t Wg ^ Wg ^ T 


The momentum effects are attributed to the fan pressure-ratio changes and pres- 
sure losses and were found to be related to the total pressure dynamics accord- 
ing to 


3p 3 m 

— = DbMp — 
3t 3t 


where b is the pressure-loss coefficient and D is a dimensionality constant. 
If we insert the liquid transport delay, the gaseous transport delay, and the 
momentum effects, equation (A29) becomes 


9p P • -TrS P • -Tr-S P • , 

— = — mre " - — mce ® + - T + DbMp — 
3t Wg Wg T 3t 


(A30) 


7pM' 


where b = 0.197 1 - I, This is equation (4) in the text. Of course, the 

T 

valve mass flow rate equations are the same, and upon substitution for T and 
arrangement the third matrix equation of figure 3 is derived. 
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Figure 1 Schematic of Langley 0.3-Meter Transonic Cryogenic Tunnel. 
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Figure 2.- Photograph of Langley 0.3-Meter Transonic Cryogenic Tunnel 





TUNNEL PROCESS MODEL 



pM 

0(1 + 0.2M^)^ 


1 +t s 
m 





e 

K.,VT(l+t SI 
M p 


p(B + a) 

i' ' '"■'l 

Ll_ 

(/> 

_/ 

w 

c 

p M 


, DbpMs '^g'’ 

P 

T0 


M 

7 3 

T6»(l + 0.2M ) 





Figure 3 


Lumped-parameter multivariable model of 0.3-m TCT 
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Figure 6 



(a) Temperature dynamics. 


diagrams of model dynamics and simulation control 
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(b) Pressure, fan, and Mach number dynamics. 
Figure 6.- Continued. 
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(c) Simulator analog block diagram. 


Figure 6.- Continued. 
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(d) Closed-loop controllers for temperature and pressure. 
Figure 6.- Concluded. 
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Figure 8.- Digital control 
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Figure 9.- Control panel for hybrid-computer simulation 
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Figure 10.- Simulator instrumentation. 
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Figure 12.- Simulator and tunnel transient response data at T = 100 K, 

M = 0.6, and p = 1.57 atm. 
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Figure 13.- Simulator and tunnel transient response data at T = 100 K 

M = 0.9, and p = 1.57 atm. 
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Figure 15,- Simulator and tunnel transient response data at T = 100 K 

M = 0.6, and p = 3.00 atm. 



TOTAL PRESSURE. Atm TEST-SECT I ON MACH NO. TOTAL TEMPERATURE , 





TOTAL PRESSURE, Atm TEST-SECTION MACH NO. TOTAL TEMPERATURE , 





TIME. Sec 

Figure 17.- Simulator and tunnel transient response data at T = 100 K, 

M = 0.3, and p = 5.00 atm. 
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Figure 19,- Simulator and tunnel transient response data at 

M = 0.9, and p ® 4.93 atm. 
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Figure 20.- Simulator and tunnel transient response data at T = 200 K, 

M = 0.3, and p = 1.57 atm. 
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Figure 25.- Simulator and tunnel transient response data at T = 200 K, 

M = 0.9, and p = 3.00 atft. 
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Figure 26.- Simulator and tunnel transient response data at T = 200 K, 

M = 0.3, and p = 5.00 atm. 
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Figure 27.- Simulator and tunnel transient response data at T = 200 K, 

M = 0.6, and p = 5.00 atm. 
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Figure 28.- Simulator and tunnel transient response data at T = 200 K 

M = 0.9, and p = 5.00 atm. 




SIMUUTOR 



TUNNEL‘S 


-100 rpm to N 





TIME, Sec 


Figure 30. - Simulator and tunnel transier 
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Figure 33.- Simulator and tunnel t: 

M = 0.6, and 
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Figure 35.- Simulator and tunnel transient response data at T = 275 K, 

M = 0.3, and p = 5.00 atm. 
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Figure 37.- Simulator and tunnel transient response data at T = 275 K, 

M = 0.8, and p = 5.00 atm. 
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Figure 38.- Cool-down profile for tunnel and simulator. 
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Figure 40.- Comparison of tunnel and simulator warra-i 


•nI 


WARM-UP CONDITIONS : 



SIMULATOR 

Mach no. held constant at 0,6 
Total pressure held constant at 2 atm 

TUNNEL 

Mach no. held at approx. 0.6 
Total pressure held at approx. 2 atm 


profiles. LN2 valve closed at 0. 



PRESSURE 
SET POINT - 


GAIN SCHEDULE 
AND 

CONTROLLER 


PRESSURE LOOP 


Figure 41 Pressure- 
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schematic for 0.3-m TCT control. 
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Figure 42. 
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ture-loop schematic for 0,3-m TCT control. 











(a) Low tunnel temperature. 

Figure 43.- Simulator temperature and pressure closed-loop responses 
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(b) High tunnel temperature. 
Figure 43.- Concluded. 
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